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Abstract 

We study the possible stationary persistent supercurrents flowing on a cylindrical sample 
supporting a two-dimensional charged fluid. The internal dynamics of the fluid is obtained 
by means of an effective theory in which the fluid self-interacts through a U(l) gauge field. 
We find that the presence of persistent supercurrents depends on what kind of gauge field 
it is. In particular the current is zero if it is a Maxwell gauge field, and it is maximal if 
it is a Chern- Simons gauge field. There is an intermediate behaviour in presence of both 
Maxwell and Chern-Simons term. Therefore it appears that persistent supercurrents are 
possible only if the fluid is chiral. 
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1. Introduction and summary, 



In this paper we study the possibility that stationary supercurrents can flow, like in a 
superconductor, in a planar non-relativistic fluid which looks incompressible in the sense 
that all the excitations have an energy gap and thus there are no gapless compressional 
modes. 

This fluid is described by the most general effective action in two space and one time 
dimensions including U(l) gauge fields. As it is well known (references [1] [2]), in 2 + 1 
dimensions the generic action with at most two derivatives for the gauge fields includes 
both a Maxwell- like term and a parity and time reversal violating Chern-Simons term. A 
gauge theory of that kind has been introduced long ago and extensively studied also for 
possible relations with Chern-Simons theory describing fractional statistics particles and 
rather speculative applications to high T c superconductors (see reference [3] for a general 
introduction) . 

Here we study an effective action where the U(l) gauge field is coupled to a non- 
relativistic matter field. This effective theory was studied for the pure Chern-Simons case 
in reference [4]. The possible vortex excitations of this system for the Chern-Simons case 
and for the Maxwell case were discussed in reference [5], and later also in [6]. (An effective 
theory formally similar to the one we investigate here, in the Chern-Simons limit, was 
used to describe the Fractional Quantum Hall Effect in references [7] , [8] , [9] . Notice that 
the physical system we consider is completely different from the quantum Hall setting, in 
particular there is no uniform magnetic field, provided by a given external device). 

One can easily derive the spectrum (see section 2) and find that the possible excitations 
have an energy gap. Thus, this fluid does not appear similar to a standard superconductor 
(nor to the standard anyonic superconductors [10], [11], [12], [13], [14]), which possesses a 
gapless excitation that plays the role of the Higgs field when coupled to the electromagnetic 
field. Nevertheless, we specifically consider the possibility of stationary currents flowing 
around in an idealized cylindrical sample, as the one of figure la, much like a supercurrent 
in a superconducting solenoid. We consider the cylinder made of a pile of many two 
dimensional annuli, on which the fluid's current flows around. We also consider the fluid 
to be charged and coupled to ordinary three-dimensional electromagnetism, and we study 
in particular the magnetic field produced by the possible supercurrents. 

We find in section 3 that the possibility of supercurrents depends very much on the 
relative weight of the Maxwell term as compared to the Chern-Simons term, in the two- 
dimensional (gauge) dynamics. An important role is played by the electrostatic force which 
arises when there is a charge fluctuation (we assume that our system, fluid plus background, 
is overall electrically neutral. However we assume that the background is rigid, and thus 
a fluctuation in the density of the charged fluid gives rise to an electrostatic field) . In the 
limit where the Maxwell term dominates, no currents are possible in the bulk of the sample 
- in this case the currents are concentrated at the edges and are small. Thus in this limit 
our sample does not work like a superconducting solenoid. 

On the contrary, in the limit where the Chern-Simons term dominates, the current can 
flow freely like in a superconductor and our system behaves like a superconducting solenoid 
(of the same thickness, of course), the resulting magnetic field and flux being the same as 
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the ones due to a super current for an ordinary superconductor. In the intermediate cases 
we find intermediate behaviours. 

Since there is a gap in the spectrum, the flow of the current is protected against 
external disturbances, in the sense that as far as the energy carried by the external sources 
is less than the gap, no interaction is possible. 

Thus, for a sizable Chern-Simons term, we find that, despite being non-standard, our 
fluid has some properties similar to standard superconductors. Actually we studied [4] 
already other properties of the non-relativistic fluid coupled to a three-dimensional U(l) 
gauge field dominated by the Chern-Simons term and we found interesting results showing 
similarities with superconductors. In particular, we studied the possible penetration of an 
external magnetic field (Meissner effect) in the bulk of an idealized sample of many two- 
dimensional layers where the fluid lies. We found an intermediate behaviour, namely if 
the magnetic field is orthogonal to the layers (that is, in a configuration similar to the one 
studied here) the system behaves in the same way as a type II superconductor, whereas if 
the the field is parallel to the layers it can penetrate much more easily, going to zero inside 
with an inverse power law, rather than with an exponential law like in the standard case. 

Thus, putting together the previous results and the ones presented in this paper, 
we can conclude that a chiral fluid (chirality being related to the parity breaking Chern- 
Simons term) shows interesting properties which could qualify it, in some respect, as a 
kind of non-standard superconductor. 

2. The effective lagrangian. 

Here we study a non-relativistic theory in 2 + 1 dimensions described by the following 
lagrangian: 

1 2 1 2 

-aex^Axd^ + ^D 2 cj> + % - (V - 0^) + (H 2 - Po) -V(\(f>\) . (2.1) 

Here <f>(x, t) is a complex field which plays the role of order parameter and is related to 
the density by the relation: 

p(x,t) = \cf>(x,t)\' 2 . (2.2) 

The covariant derivative is Di = di — iAi. Notice that A Q is coupled to dp = p — p , 
and that p is the average density representing the neutralizing background necessary for 
the consistency of the theory. Since for conservation of the total number of particles the 
integral over all the surface of dp is zero, we can write it as the divergence of some vector 
field u{x) which we choose to be irrotational V A u = 0: 

5p = V-u (2.3) 

Equation (2.1) is the most general non-relativistic gauge invariant lagrangian with minimal 
coupling to a vector potential, once the maximum number of two derivatives is required. 
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First note that beside the Maxwell term there is also a Chern-Simons term with coupling 
constant a. Notice also that the coupling constants for the electric and the magnetic part 
of the Maxwell term are different. They are equal only if the theory is Lorentz invariant, 
that is considering a relativistic theory (here we draw the reader's attention to the fact 
that in 2 + 1 dimensions these coupling constants have the dimension of a length and a is 
dimensionless). In the following, for definiteness, we take V(|0|) = ^A (|</>| 2 — p ) . 
Notice in particular that varying (2.1) with respect to A one gets (in the gauge 
V -1=0): 

g B d (doAi - diA ) + g B eijdjB = J t - 2ae ij (d j A - d Aj) (2.4) 
g E A Aq = 5p — 2aB . (2.5) 

We have defined: 

8p = p-p B = e lJ 8 l A J Ji = (0^0 - <9*0V - 2iAi<^<f>) . (2.6) 

Thus for a = we recover the two-dimensional Maxwell equations (though with two differ- 
ent couplings). For g E =0, on the other hand, from equation (2.5) we get Sp — 2aB which 
is the usual Chern-Simons constraint relating the field strength to the matter density. 
Later on we will also consider the fluid to be electrically charged and thus coupled to the 
true, three-dimensional electromagnetic field. 



2.1. The spectrum. 

Now we will study the energy spectrum of the small field perturbations of our Maxwell- 
Chern-Simons theory. 
We take the following parameterization: 



= fe ie • 



(2.7) 



and rewrite lagrangian (2.1) keeping only quadratic terms in the fields. With this choice, 
and in the gauge V-A = 0, equation (2.1) becomes (in this gauge one can put Ai = €ijdjifj 
and see that J d 2 x AidoAj =0): 



C=-g E 



(d Q A t ) 2 + («Mo) 2 ] - \g B (d t A J ) 2 + 



-2ae ij Aod i A j + 



1 

2m 



1 

4p 



(d t 5p) 2 -Po (diOy-poA 2 



+edo(5p)+A Sp--X(5p) 2 . (2.8) 



From this lagrangian we get the following equations: 

5C 



SA 



—g E A A — 2aeijdiAj + bp = 



(2.9) 
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— = -g E dlAi + g B A Ai - 2ae ij d j A - —Ai = (2.10) 

S Jl = P± A e + d (8p) = (2.11) 
ov m 

a r i 

ASp-d 9 + A -X8p = 0. (2.12) 



6(5 p) Amp Q 
Multiplying (2.10) by e ki d k we get: 

-g E d$B + g B AB + 2aAA -^B = 0. (2.13) 

m 

Now multiplying equation (2.11) by <9 and equation (2.12) by A we get 



and 



Ae = --d 2 (5p) , (2.14) 

Po 

— A 2 Sp- A6 + AA -XA6p = . (2.15) 
4mp 

We can now eliminate 9 and A from equations (2.13) and (2.15) by using (2.9) and (2.14). 
By taking B = (Boe l ( £t ~P' x ^ + c.c.) and bp = (5p e l ( £t ~P' x ^ + c.c.) we thus obtain: 

£ 2 B - ( + — ) B o + ^Po - -(pl+P 2 y )Bo = (2.16) 
V 9l rng E J g\ g E v 

£ 2 S Po - ^5p + 2 -^B - ^(P 2 x +P 2 y )5p - -^(pI +p 2 y ) 2 5p () = . (2.17) 
mg E mg E m y Am z y 

Therefore the values of the energy for p = are the eigenvalues of the above system: 



g, & + _£!_ ± & Jl + !HZ. (2.18) 

9 E rng E g% \ ma 2 

Notice that Sf 2 are both positive. If a^O, that is if there is no Chern-Simons term, the 
two eigenvalues are equal: 

£i=S2 = tP^. (2.19) 



mg E 

If, on the other hand, we take the opposite limit g E -^0 then only one of the eigenvalues 
remains finite whereas the other one diverges: 



c 2 I«I 

t\ = > oo 



£■ 
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g e 

Po 

2m\a\ 



(2.20) 



That is, in general there are two propagating degrees of freedom each with its energy gap. 
In pure Maxwell theory the energy gaps are degenerate. This degeneracy is resolved by 
adding the Chern-Simons term. If there is only the Chern-Simons term the theory possesses 
only one degree of freedom with an energy gap, the other one being frozen. This value of 
the energy gap for the pure Chern-Simons case has already been derived in reference [4]. 



3. Stationary currents 



In this section we study the possibility of stationary currents. We would like to study 
a configuration that, although idealized, is near to a setting relevant for a realistic device 
which could carry superconducting currents. Therefore we take the surface where our two- 
dimensional fluid lies to be an annulus of radii r<i >r\ such that r2 — r\ =L. We consider 
a pile of many of these annuli separated by a spacing d, see figure la. We assume that 
L/Vi <C 1, but still L is macroscopic, that is L/d^> 1 where d is any length of the order of few 
atomic distances such as the inter-annuli spacing. Furthermore, the total length L z of the 
cylinder is supposed to be very large as compared to the radii r\^- We study the possible 
configuration where the fluid's current flows around in each annulus, see figure lb, and it 
is uniform with respect to z, that is it is the same in each annulus. Thus our cylindrical 
configuration will act as a solenoid. Later we will compute the three-dimensional true 
magnetic field inside the solenoid, due to the current. 



From the two-dimensional point of view the annulus is equivalent to a strip which is finite 
in the a; direction, —L/2<x<L/2, and periodic in the y direction, 0<y<2nRM, see figure 
lc (with Rm = ri \ r2 , remember that in this approximation Lj Rm <C 1). We consider a 
possible current induced by a quantized phase for the matter field, as for a vortex: 



where p= n integer, is fixed. The current flows around the cylinder as in figure lb. 



Notice that we consider a configuration with cylindrical symmetry, therefore we take (in 



Figure 1 



(j) = \4>\e ipy 



(3.1) 



the gauge V • A = 0) : 



u = 



A= (A x , A y ) = (0,A(x 
(u(x),0j 5p(x) = d x u(x) 



)) B = d x A 




(3.2) 
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As we said we assume uniformity in the z direction and we consider the limit where the 
cylinder's length L z is very large as compared to Rm, so we can write the hamiltonian of 
the model as: 



+ ^-(Sp- 2aB) * j^— * (5p - 2aB) . 

2g E (-A) 



(3.3) 



Here the last term comes from integrating out the auxiliary field Aq and the symbol * 
means convolution. 

It is convenient to express this hamiltonian as the integral of a local hamiltonian density. 
To this aim we rewrite the last term of (3.3) in terms of the two-dimensional "electrostatic" 

— * 

field E determined by the equation: 

V-£ = - — (5p-2aB) (3.4) 
9e 

and imposing boundary conditions such as to eliminate possible zero modes which do not 
contribute to the hamiltonian (3.3). This amounts to requiring (in our configuration with 
cylindrical symmetry we have uniformity in the y coordinate, therefore E = (E(x), 0)): 

E(x < -L/2) = E(x > L/2) = . (3.5) 

The hamiltonian can then be rewritten as 

|- = | dxdy \^g E E 2 + \g B B 2 - ^ D 2 <j> + V (\<P\)} (3.6) 
We have to fix completely the boundary conditions. We assume that: 



A(x < -L/2) = A(x > L/2) = (3.7) 

that is, we assume A, which represents in general a propagating excitation of our fluid, 
to be zero outside the sample i.e. where there is no fluid. Notice that, by subtracting the 
irrelevant constant cq = —\\p N where iV is the fixed total particles' number N = J |</>| 2 , 
we can rewrite J V=^X J|</>| 4 which vanishes outside the sample. 

One can also say that, by definition, dp = outside the sample. Since the fluctuation of 
the particles' total number in a given domain is the flux of u through its boundary, and 
the particles' total number is fixed, the boundary conditions for u are 

u(x < -L/2) = u(x > L/2) = (3.8) 



Thus, keeping into account (3.5), we see that the hamiltonian density vanishes outside 
the sample, consistently with the fact that it describes the dynamics of the fluid which is 
confined in the strip. 
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We can check that the boundary condition (3.5) is the one which correctly ensures that 
the hamiltonian (3.6) implements the Chern-Simons dynamics in the limit g E ^0. In fact, 
the solution of equation (3.4), with the boundary conditions (3.5), (3.7) and (3.8) gives: 



X 



E(x) = -— [ dx'Spfx') + —A(x) = -— \u(x) - 2aA(x)] . (3.9) 
9e J 9e 9e 

_ L_ 

2 

Thus we see that, since finite energy requires g E E 2 to be bounded, one recovers in the 
limit g E ^0: 

A(x) = (3.10) 

which is indeed the solution of the Chern-Simons constraint B — j^Sp with the boundary 
conditions of (3.7) and (3.8). 

To this hamiltonian we now add one extra term which plays an essential role in what 
follows. This term represents the true, three-dimensional electrostatic interaction between 
the fluctuations of the charged matter. Of course, the whole system, fluid plus background, 
is overall electrically neutral. But if the background is rigid, as we assume here, a local 
density fluctuation of the charged fluid gives rise to electrostatic forces. This yields a real 
electric field which obeys to the three-dimensional Maxwell equation: 

V • E em = e5p( 3) , (3.11) 

where 5p^ = 5p/d is the three dimensional density and e is the electric charge, (for 
other effects of this electrostatic contribution we refer to [4]). We can then compute this 
electrostatic contribution. The three-dimensional Maxwell equation can be written as (in 
our uniform in y and in z configuration): 

d x E em (x) = e -5p= 6 -d x u . (3.12) 



This equation can be integrated yielding: 



E em = -u (3.13) 

(Jj 

Therefore the electrostatic contribution to the hamiltonian is 

d -(E-f = y/ (3-14) 

Thus, after taking everything into account, the hamiltonian of the system can be written 

as: 



^— = fdx{ 
RmL z J { 



1 .(„ _ 2aA) 2 + \g B {d x Af + J-(p - Afp + J-(«9 X H) 2 + 



2tyRmL z J 1 2g E 2 2m ' 2m 



+2*(0*«) +Yd u ) (3 - 15) 

Before going on it is maybe worth to stress again the difference between E(x) and E em (x). 
The first one is the two-dimensional "electric field" coming from the internal dynamics 
of the fluid. The second one is the actual, real, three-dimensional electric field due to 
electrostatic interaction of the charged matter. 
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3.1. A heuristic analysis. 



We perform a heuristic analysis by considering a situation where the matter is uni- 
formly distributed and there is no "magnetic field" , that is: 

dp = B = d x A = . (3.16) 

This, in particular, implies p = p and d x \<j)\ =0. In this case the variation of the hamiltonian 
(3.15) gives the equations 

5A\2irR M L z ) g E y > m y P > 

8 f H \ 1 , „ N e 2 { } 

— — - — — = — (u - 2aA) + — u = . 
ou \2ixR M L z ) g E d 

The solution of these equations is (in this heuristic analysis we forget the boundary con- 
ditions): 

Po 

A = m 

Po , 4a 2 e 2 

m d + e 2 g E ( 3 - 18 ) 
2ap d 

fe 

This corresponds to the current density: 

2„2 



p (d + e 2 g E ) + 4ma 2 e 2 ^ 



Aa 2 e 



J=^ d + f 9 2 E 2 P. (3.19) 
m p Q 4a z e z 



m d + e 2 g E 

Notice that if e 2 = 0, that is if there is no electrostatic effect, the density current is zero. 
Notice also that in the Maxwell limit (a^O or g E ^oo) the current density is still zero. 
Therefore the current density is different from zero only if there is some Chern-Simons 
amount in the fluid's dynamics. In particular for the case of pure Chern-Simons, with 
ct^oo, we recover the maximum value for the current: 

Po 4ma 2 e 2 p . , 

J — > — > —p . (3.20) 

gs^o m p d + Ama z e z o^oo m 



9 



3.2. A more detailed analysis. 



The assumptions of the previous section are physically too restrictive because the 
actual physical situation can have dp and B different from zero. Nevertheless in this section 
we show that a more accurate computation leads substantially to the same results, apart 
from a correction which is different from zero only in regions very close to the sample's 
borders. 

In order to keep the discussion reasonably transparent, we make the simplificatory assump- 
tions that p ^>\Sp\ and that l/m<gA. With these assumptions we can write: 



p{p~Af ~ Po (p-Af 



i 



8mp 



(djp) 2 ~ 



(3.21) 



We stress that our results would hold also in general. In fact, the equilibrium configuration 
of the fluid, far from the borders, will be the one of section 3.1. What happens in details 
near to the borders will depend on the details of the hamiltonian, and we have chosen to 
present the simplified discussion based on (3.21). 
In this way the field equations now read: 



_5_ 



H 

2itRmL z 
H 



= -d£A + 



2a 



(2aA-u) + — (A-p) = 



-d 2 x u+ , 



9e9b 
1 



(3.22) 



(u - 2aA) + —u 



5u \2tvRmL 

These equations can be rewritten in matricial form as follows: 

(-10* + M)tf = $ 
where 1 is the unit two by two matrix, and 



(3.23) 



M 



( 



4a z 



+ 



Po 



2a 



\ 



9e9b mg B 
2a 



\ 



9eX 



9e9b 

_L i_ 

9e A Xd 



/ 



Po 

g B m' 





-p 



Note that detM^O. 

The general solution of equation (3.23) can be written as 



where 



A* 



(3.24) 



(3.25) 



(3.26) 



is one particular solution and A* and u* are given by equation (3.18). \l/o is the general so- 
lution of the homogeneous equation (— Id 2 + M)\I> = 0. The boundary conditions discussed 
in section 3 are L/2) = \I>(L/2) = 0. Therefore we can write: 



A(x) 
u(x) 



cosh( / Uix) / A\ \ cosh( / u 2 a:) 
cosh(^iL/2) \ ux J cosh(^ 2 -^/2) 



A 2 
u 2 



+ 



A* 
it* 



(3.27) 
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where 

A»i,2 = g ( a + d ± ^(a - d ) 2 + 4bc ) ( 3 - 28 ) 
are the eigenvalues of the matrix M, and 

m,2 = //l ' 2 b ~ a Ai, 2 ■ (3.29) 

One can verify that //f 2 are rea l positive. Let us mention their values in some limiting 
cases. For the case of pure Maxwell, i.e. a = 0, the fields A(x) and u{x) decouple and the 
two eigenvalues are 

n1 = — — (relevant for A) ^ = ~^ (relevant for u) . (3.30) 

In the opposite limit g E ~^0 the relevant field configuration is u = 2aA, and this is reflected 
in the eigenvalues which turn out to be 

Po 4e 2 a 2 
— + 



Imposing the boundary conditions, we determine A\^ by: 



1 1 



b b 




ri-a ^-a ) ( 7 )=-( 7 ) • (3.32) 



Notice that the two eigenvalues can be equal only for — = and -p— = — + -fj in which 

° " " gE rngg ggA Ad 

case M is proportional to the identity, and in (3.27) we can write A\ = A 2 and u\ = U2 
with no relation between A\ and u\. The boundary conditions give in this case A\ = — A* 
and u\ = —u*. 

Since n\ 2 are expressed in terms of the microscopic parameters of the hamiltonian, we 
have L-/ii,2 ^> 1. Thus, A(x) = A* and it (a;) = u* apart from very near to the edges, that is 
for |x| ~L/2, and therefore we recover the solution found in the previous heuristic analysis. 



3.3. Magnetic field. 



In this subsection we compute the actual magnetic field present inside the cylinder 
(see figure la) due to the supercurrent computed above. 

It is clear that we have to take now a three-dimensional view-point. Our three-dimensional 
system is now the cylinder embedded in the three-dimensional space. Therefore we will use 
the cylindrical coordinates (r, (p, z) as in figure la. We remind that we consider L<^ri j2 so 
that the cylinder can be thought of zero thickness and we take its radius to be Rm = ri t r2 . 
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The three-dimensional density of the fluid can be written in this approximation as p^ = 
fL-8(r-R M ). 

Furthermore we suppose L z ^> Rm so that it is possible to keep far from the cylinder's 
edges and therefore to disregard the edge effects. 

With these definitions the three-dimensional current density flowing all round the cylinder 
is: 

J v = ^[P~ e A 7( r ) ~ A ] 6 ( r ~ R m) ■ (3-33) 

Here A < ^ n (r) is the electromagnetic vector potential describing the magnetic field produced 
inside the cylinder. Since we limit ourselves to the region far from the edges of the cylinder, 
we can assume A e ™ independent of z. 

A is the (two-dimensional) vector potential describing the fluid's dynamics on the strip 
considered in the previous section. 

We can reproduce all the computation of the previous section taking into account for A^ 
simply with the substitution p-^p — eA e ™. Therefore we now have, instead of (3.18): 



Po 

m 



2„2 



Po 4a 2 e 



■(p-eA™) (3.34) 



m d + e 2 g E 

We have neglected the deformation of A very near to the annuli's edges, as discussed in the 
previous section. Now we use this expression for A to solve Maxwell's equation V AB em = J 
inside the cylinder (see reference [15]). The solution of this equation is: 



ept™ 1 



m e p Km 

2m 

where: 



-e(R M - r) + ^6(r - R M ) 



(3.35) 



A 2 2 

4e a 



PT= d + e " g 2 E 2 -Po§- (3.36) 
p 4e a a 

m d + e 2 g E 

Here 6 is the step function. This yields the current density: 

J, = Am E S(r - Rm) , (3.37) 

2m~ 

We can also compute the magnetic field inside the cylinder using the cylindrical coordinate 
relation B e z m = ±£ {rA e ™) : 

B em = P q{Rm _ r) (3 38) 

m e 2 p^ m R M 
2m 
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This result is the same as the one of a solenoid in which the current density (3.37) flows. 
We can also compute the flux of the magnetic field in the cylinder: 



em i 

HB em ) = nRlBr = nR 2 M ^ p . (3.39) 

2m 

Remembering that p = we see that for epl m R 2 M — >oo we get 

$(5—) _^ . (3.40) 

e 

That is, the total flux is quantized in the above limit, as expected for a vortex-like current. 
We stress that, as we said, the result for the magnetic field and flux is exactly the same 
as it would have been obtained for an ordinary superconductor flowing on a cylindrical 
surafce, p^ m being its surface density. (It is a general fact for any superconducting current 
flowing on a cylindrical surface, that the flux quantization is strictly speaking obtained 
just in the above limit). 

Note that pl m vanishes in the Maxwell limit g E — > oo (or a = 0) and p e ™ = p j[ for 
a — > oo. Therefore we see again that, in order that our idealized device could work as 
a superconducting solenoid, there cannot be only a Maxwell dynamics, but the chiral 
Chern-Simons term must play an essential role. 
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Figure caption 

Figure 1. 

a) The cylinder pile of many annuli 

b) A single annulus 

c) The annulus as a periodic strip 
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